Drell-Yan plus missing energy as a signal for extra dimensions 
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We explore the search sensitivity for signals of large extra dimensions at hadron collid- 
ers via the Drell-Yan process pp £~^£~ + fJx X {£ = e, /x) where the missing transverse 
energy is the result of escaping Kaluza-Klein gravitons. We find that one is able to place 
exclusion limits on the gravity scale up to 560 GeV at the Fermilab Tevatron, and to 4.0 
(3.3) TeV at the CERN LHC, for n = 3 (4) extra dimensions. 
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I. INTRODUCTION 

on" 

CTN ' The idea of the existence of extra spatial dimensions has fascinated physicists for nearly a century A 
quantum theory of gravity seems to be consistent only via extension to extra dimensions Q| . If compactification 
of the extra dimensions occurs near the Grand Unification scale or the Planck scale 10^^ — 10^* GeV, then the 
effects of quantum gravity would be accessible only at very high energy scales, beyond the reach of any collider 
^ ' experiment. In this case, one would have to understand the mechanism of weak scale stabilization (at about 
100 GeV) against radiative corrections, the so-called hierarchy problem. Recently, a radical scenario has been 

■ advocated wherein quantum gravity may become significant at a much lower energy scale (Mg), as low as 
[ ©(TeV). The apparent large Planck scale {Mpi ^ 10^^ GeV) is then attributed to the more rapid l/r"+^ decrease 
. of the gravitational force with distance in n extra space dimensions. In terms of the large compactification size, 

CN ' R ^ ^/Ms, of the extra dimensions this leads toQ 
> ' 

Cn . M^i ~ R'^M'^+^ . (1) 

' Such a scenario alleviates the hierarchy problem by restating it as a geometrical one, namely understanding the 
lO \ size of the compact dimensions and the scale Mg at which compactification occurs. This scenario could be realized 

■ in certain string formulations 0. Naturally, after compactification, there will be towers of Kaluza-Klein (KK) 
excitations with mass separation of 0{1/R). To avoid conflicts with the Standard Model (SM), it is assumed that 
the SM fields are stuck to a 4-dimensional hyper-surface, while only gravitons propagate in the extra dimensions. 

' If we are interested in low-scale quantum gravity effects with Ms ^ ©(TeV), the minimal scenario n = 1 has been 
Qh, ruled out because R '-^ 10^ km, which would yield effects visible in planetary motion. For n = 2, although there 
^ is no conflict with Newtonian gravity or astronomy for R ^0.1 mm, the constraint from supernova cooling has 
, put a bound on the scale up to Ms > 50 TeV M. One may even be able to push the scale up to Ms > 110 TeV 
from the spectrum for diffuse gamma radiation^S| . For n > 2, there is no experimental or observational conflict 
^ with the theory. 

Most interestingly, there will be significant consequences for low-energy phenomenology with this scenario. 
Although the coupling for an individual KK excitation is gravitationally suppressed, the cumulative effect from 
^ , the tower of states is suppressed only by l/M^"*"^. Generically, there are two classes of collider signals induced by 
the KK gravitons. First, real KK gravitons can be emitted off SM fields. Second, virtual KK gravitons may be 
exchanged between external SM fields. There have been many studies proposing searches for visible signatures 
of extra dimensions at colliders P,|6|,p^-p^ . In this Letter, we study another process to explore the sensitivity 
to probe low-scale quantum gravity effects, Drell-Yan charged lepton pair production, which has previously been 
used as a powerful test of the Standard Model and has provided severe constraints on new physics. For the process 
of current interest, we look for the clean signal 

PP ^ e+r + pT X, (2) 

where the missing transverse energy is due to an escaping KK graviton. We calculate the signal rate and those 
of the leading SM backgrounds, and study the sensitivity to probe quantum gravity effects both at the Tevatron 
and at the LHC. 



X 



^The precise relation depends on the convention for R and AIs- It is taken to be Gjy^ = 47r i?"Mg^^ in and 
to be Gjf^ = 8TrR"-Mg^^ in where Gjj^ = M^^ is the Newton's constant. For the convenience of calculations 
for physical quantities, we take the normahzation to be [| G^^ = i?"M^+V(47r)"/2r(n/2). 
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II. CALCULATIONAL TOOLS 



The signal subprocess is 

qq -> i+r + Gkk (3) 

where £ = e, /j,, and the graviton Gkk escapes the detector, resulting in missing energy. The signal can be 
described by fourteen tree level Feynman diagrams, seven each for Z and 7 exchange, where a graviton is attached 
to each SM field and also to each SM vertex. Thus, off-shell effects of the Z boson are fully included. This would 
be extremely tedious to calculate by hand, or even with the aid of a trace-based Feynman graph program, 
but we can make this straightforward by using the helicity amplitude method: summing the numerical values 
of individual Feynman graph amplitudes for a set of fixed external helicities and momenta, then squaring the 
summed amplitude and integrating over all possible helicities and phase space numerically. 

To do this, we constructed three new HELAS vertex routines, for the graviton- fermion-fermion, graviton- 
gauge boson-gauge boson and graviton- fermion-fermion-gauge boson vertices ||^ . Coding of the matrix elements 
can then easily be done by hand. We have verified current conservation for the full matrix elements both 
analytically and numerically. The summation over KK states with different mass is carried out numerically based 
on the weight function in 

As we are including photon interference effects, we must allow the Z to be off-shell, and thus must include 
finite-width effects for the Z propagators. This presents a problem for the calculation, as the graviton amplitudes 
are not gauge invariant for a Z propagator including an imaginary piece. While we believe we can provide a 
formal prescription for properly including such a term, for the moment we must rely on an approximation that is 
known to be extremely reliable jist : we do not include a finite width in the matrix elements, but instead multiply 
the summed-squared amplitude by an overall factor 



{§ - M|)2 + (MzTz)^ {m^ - M|)2 + {MzTzy 



(4) 



where s is the parton cm. energy squared and niu the invariant mass of the lepton pair. This factor removes the 
zero- width propagators from the matrix elements and inserts Breit-Wigner resonances for the Z boson. For phase 
space regions where the principal contribution comes from graphs where the incoming quark pair or outgoing 
lepton pair have invariant masses far from the Z mass, this factor is essentially unity. When either pair is at or 
near the Z mass, it approximates the full cross section by the correct resonant contributions. 



III. RESULTS AND DISCUSSION 

For the numerical evaluations, we use CTEQ4L parton distribution functions pq] and the EW parameters 
mz = 91.19 GeV, m* = 175.0 GeV, sin^ Ow = 0.2315, and Gf = 1.16639 x 10^ GeV^^ choose the 

factorization scale fJ-f = Et of the lepton system. We impose basic acceptance cuts for event identification, based 
on detector capability. These are 

PT{i) > 15 GeV, \r]e\ < 2.0 for Tevatron, (5) 
PT{e) > 20 GeV, \r]e\ < 2.5 for LHC, (6) 

where prii) and rji are the transverse momentum and the pseudo-rapidity of a charged lepton. 

A defining feature of the signal is the missing transverse momentum due to the escaping massive graviton. Also, 
most signal events are produced in association with a Z boson. The irreducible SM background to this signal 
is from events, which are dominated by Z^*'>Z and 7*Z production, which we call "Drell-Yan-t-z/i/" . We 

want to reduce the photon continuum from ^* Z events where the virtual photon produces a relatively low mass 
lepton pair, and therefore require 

> 20 GeV, m« > 10 GeV . (7) 

For the Tevatron at = 1.8 TeV (Run I) with these cuts, this SM background is 9.4 fb; at = 2.0 TeV (Run 
II) it rises to 11 fb. The slightly higher \/j and greatly increased luminosity expected for Run II allow us to 
impose an additional cut on the missing transverse energy, 

> 100 GeV , (8) 
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FIG. 1. Normalized distributions of the lepton pair at the LHC; (a) azimuthal opening angle 4>u; (b) ARu. 
The n = 3 graviton signal is shown with the solid curves, and the SM background with the dashed curves. 



reducing the SM background to only 1.25 fb. 

The much higher Vi available at the LHC allows for the emission of much heavier gravitons and for a significant 
recoil of the dilepton system. This is reflected in the normalized dilepton angular distribution of Figure |l|, for 
an n = 3 signal and the DY+vv background: (a) azimuthal opening angle and (b) separation ARu — 
+ ivti — %2)^- These plots show distinct differences for the signal and background. The final state leptons 
of the signal are preferentially emitted in the same direction, close to each other, while in the background the 
leptons tend to be more back-to-back. We can thus impose further cuts on the leptons, requiring 

<^«<90°, ARa<1.2. (9) 

This heavy graviton emission will also result in the signal exhibiting a much harder fiJx spectrum, as shown in 
Fig. H We therefore suggest an even higher cut for the LHC, 

^T>150GeV. (10) 

After the cuts of Eqs. (|,|J|,|l^), the SM background at the LHC is 2.8 fb. 

A feature of any effective theory is that it is a low-energy approximation and can not be trusted at large 
energies, comparable to the defining scale of the theory. Applying this rule to the present case, our calculation 
requires some additional care for n > 2, as a non-trivial fraction of those signal events occur at center of mass 
energies \/I > Ms- A conservative estimate of the signal is obtained by discarding any events with ^/l > Ms 
More generally one could invoke form-factor damping of the high energy region. To obtain estimates of the 
string scales which can be probed experimentally we employ an iterative procedure, starting with a seed value for 
Ms and throwing away all events with \/I > A/5. The resulting signal cross section is used to recompute a new 
exclusion limit for Ms, assuming a simple scaling behavior of the signal cross section, a signal ~ l/Mg"*"^. Using 
the new exclusion limit as the upper bound on the allowable -s/l range, one obtains a stable exclusion limit in 
only 2-3 iterations. However, for n > 4 most of the signal events populate this unphysical region, and the overlap 
of the validity range of the effective theory with the range which would provide for a visible signal shrinks to 
zero. In this situation, one is unable to obtain any reliable estimate for a probe to Ms via signals. Instead, 
one would expect to first observe string excitations of particles states as a signature for new physics JT^ . 
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FIG. 2. Normalized missing transverse momentum distributions of the n = 3 graviton signal (solid) and SM 
background (dashed) for the LHC. 



n=2 n=3 n=4 



Tevatron Run I <e 


0.2 fb"' 


0.9 






Tevatron Run II ( 


a 10 fb"^ 


1.5 


0.56 




LHC @ 100 fb"^ 




5.3 


4.0 


3.3 



TABLE I. String scale Ms 95% CL exclusion limits (TeV) at the Tevatron and LHC for n ^ 2,3,4 extra 



dimensions. For the Tevatron Run I value, we have imposed the cuts of Eqs. (|5|, 7 
Eq. (^) for the Tevatron Run II, while for the LHC the cuts are given by Eqs. (3 
account the a/J < Ms requirement discussed in the text. 



and the additional cut of 
flCl). The limits take into 



In Table we summarize the sensitivity reach to the scale Ms for n = 2 — 4 via the process of Eq. (|^) at the 
Tevatron and the LHC. We see that this mode is a quite promising search channel. Although the Tevatron can 
impose some bounds, it is much more impressive to search for the signal at the LHC, its reach being several TeV 
for n = 3, 4 extra dimensions. 

A comparison with other studies at hadron colliders is in order. While the monojet+^T signal from 
processes like qq ^ g Gkk has the largest rate, it also has much more severe QCD backgrounds, mostly due to 
the mismeasuremcnt of jets in the forward regions of the detector. Ref. fl^] obtained results comparable to ours. 
For instance,^ at the LHC a 95% CL limit is expected for Ms = 6.4 (33) TeV with n = 2 (4). On the other 
hand, Ref. [g| reached a more impressive conclusion, claiming a 5a discovery for AIs = 14 (6.0) TeV with n = 2 
(4) extra dimensions^. Alternatively, the virtual Gkk contribution to the DY process of qq can be also 

significant. It was found that a 95% CL limit can be reached a scale A ~ 1 (6) TeV at the Tevatron and the 
LHC pr] , with a little dependence on n. However, one would have to introduce an additional assumption for a 
cutoff scale Aabove which the virtual KK tower is truncated, making a direct comparison of Ms exclusion limits 
from external Gkk production versus virtual exchange difhcult. 



^We have converted the scale Ms to our normalization convention. 
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IV. CONCLUSIONS 



We have explored the search sensitivity for signals of large extra dimensions at hadron colliders via the DrcU-Yan 
process pp £^£^ + X {(. = e, fj) where the missing transverse energy is the result of escaping Kaluza-Klein 
gravitons. This is a very clean channel for hadron collider physics. We find that one is able to place exclusion 
limits on the gravity scale up to 560 GeV at the Fermilab Tevatron, and up to 4.0 (3.3) TeV at the CERN LHC, 
for n = 3(4) extra dimensions. This reach is comparable to one found in previous studies of the monojet+jE'T 
process [M and for virtual contributions of Gkk to the DY process [O . 
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